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FOREWORD 


The  work  described  in  this  report  was  authorized  under  Project 
1B522301A081 ,  Chemical  Agent  Dissemination  (U) .  This  work  was  started 
in  June  1966  and  completed  in  May  1967  for  thia  reporting  period. 


Reproduction  of  thia  document  in  whole  or  in  part  is  prohibited 
except  with  permission  of  the  CO,  Edgevood  Arsenal,  ATTN:  SMUEA-RPR, 
Edgewood  Arsenal,  Maryland  21010;  however,  Defense  Documentation  Center 
is  authorized  to  reproduce  this  document  for  United  States  Government 
purposes . 


The  information  in  this  report  has  not  been  cleared  for  release 
to  the  general  public. 


DIGEST 


The  program  objective*  are  to  develop,  to  design,  and  to  £abricate 
«  complete  system  for  tha  assessment  of  aerosols  dispersed  in  rest 
chambers,  The  main  activities  reviewed  are  aerosol  mass  decay  analysis, 
chamber  mixing  and  material  balance  studies,  slsctrostatic  effects  on 
chamber  aerosols  and  their  evaluation,  bioef fectiveneie  of  seroeola, 
sampling  and  subsequent  data  reduction  methods  requited  for  the  deter¬ 
mination  of  aeroaol  parameter*  through  the  mesa  decay  procedure.  The 
design  of  an  automated  sequential  filter  sampling  system  is  described  as 
wall  as  initial  phaso*  in  the  development  of  a  beta  abaorption-lmpactor 
instrument . 

Experiment*  in  the  above  area  were  conducted  both  in  a  large  (213  m^) 
teat  chamber  with  exploaive  diaaemination  of  agent  and  simulant  materials, 
and  in  small  laboratory  chambers  with  pneumatic  aerosol  injection. 
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AEROSOL  DISSEMINATION  ASSESSMENT 


I.  INTRODUCE  ION 

Basic  program  objectives  which  have  been  outlined  in  previous  reports 
(Ref.  1-7)  are  restated  here  for  inroediate  reference. 

The  objectives  of  this  program  are  to  develop,  design,  and  fabricate 
a  complete  chamber  aerosol  assessment  system.  The  function  of  this 
system  is  to  provide,  by  the  selection  and  Integration  of  appropriate 
instruments  and  analytical  methods,  a  detailed  description  of  chamber 
aerosols  from  as  near  as  possible  to  the  instant  of  dissemination  to  time 
periods  up  to  1  hoar.  By  continuous  or  incremental  sampling,  the  state 
oi  the  aerosol  should  be  definaM.e  over  time  intervals  aa  small  as  1  second. 

Sin'r  th-  sy.  em  to  be  use.'  as  .  ma  .  tool  in  the  development 
of  chemical  device!  ,  it  should  be  adap cabl r  to  the  assessment  if  a  broad 

range  of  disseminators  and  agents  and  tc  mot  test  chambers  now  used  by,, 
groups  within  Edgewt  .«->  Arsenal  and  by  contractor  organisations .  Further¬ 
more,  extraneous  efi* cue  attributable  to  chamber  geometry,  aerosol  mix¬ 
ing  ,  wall  losses,  electrostatic  chargs ,  and  agglomeration  should  not 
interfere  with  proper  assessment  of  disseminator  functions. 

For  a  broad  capability  in  aerosol  assessment,  component  instrumen¬ 
tation  should  permit  the  following  measurements i  mass  concentration  of 
active  agent,  concentration  of  degraded  agent  and  inert  mater  tale, 
descriptive  particle-size  parameters  for  active  agents,  en<l  shape  factor 
tor  solid  particles.  In  addition  to  a  broad  rango  of  particle  size,  from 
less  than  lu  to  ICOu,  particle  dlmenalons  should  be  relstable  to  thosa 
describing  biological  effectivoness  of  aerosol  clouds. 

In  addition  to  the  requirement  that  system  components  should  define 
aerosol  cloud  properties  described  above,  it  is  desirable  that  sensor 
Outputs  from  these  Instruments  should  be  amenable  to  direct  display  and/ 
or  data  procsssing  and  storagt.  The  extent  to  which  date  are  displayed 
or  processed  and  scored  for  later  recall  will  depend  upon  vnich  instruments 
or  combinations  thereof  best  meet  the  system  objectives.  In  certain 
lnetancas,  automation  may  be  highly  impractical  on  the  basis  of  the  cost 
of  the  process  and  the  accuracy  with  which  measurements  can  be  made  with 
available  equipment.  Since  the  aerosol  assessment  system  is  to  be  fur- 
nlshtd  as  a  tested,  prototype  assembly  at  the  end  of  a  3-year  period, 
the  cnntiact  framework  precludes  extensive  exploration  of  new  or  novel 
instrume  .tal  concepts.  Therefore,  it  hes  been  recommended  that  com* 
marcially  available  equipment  or  well-advanced  experimental  devices  be 
used  ac  the  basic  system  components  and  that  these  components  oe  modified 
and  calibrated  in  accordance  with  specific  system  needs. 


II.  REVIEW  i ip  SECOND  YEAR  P '.mV. RAM 


The  techniques  and  ins  t  rira  .it  s  w  ich  /.re  to  be  in  -  luded  in  the 
chamber  assessment  system  were  discussed  i.i  "he  Sixth  oeavtcrly  Ptogress 
Report  (Ref,  6).  To  recapitulate,  the  basic  aerosol  measurement  concept 

tail  i  jo  A  \j  y  a  h  l  i.  k  C  1  CUitij'  i  iiif,  n  s  l  n  c  i-  cut  l  fji  u  »  tmv  w  <.  *  vt  *  nwnii-cu  *-  uc 

of  the  collected  mass,  Chemical  analysis  is  to  be  uued  to  discriminate 
between  active  agent  and  background  components.  Particle  siting  is 
based  upon  mass  decay  curve  inversion,  cascade  impaction,  light  scat¬ 
tering  particle  counting  as  a  function  of  time,  and  microscopic  analysis 
of  samples  on  membrane  filters  and  other  plane  burfaces .  The  latter  is 
also  etn ployed  to  distinguish  shape  factor.  Chemical  analytical  methods 
are  to  be  implemented  by  automatic  analysis.  A  variety  of  other  measure¬ 
ments  are  to  be  made  to  support  the  test  analysis  including  chamber 
temperature,  pressure,  humidity,  electrostatic  effects,  etc. 

A  number  of  activities  have  been  pursued  during  the  past  year  to 
establish  the  validity  of  the  selected  assessment  methods  and  to  modify 
the  me  the  is  for  use  in  the  system. 

Mass  decay  inversion  analysis  was  pursued  intensively  since  this 
technique  has  been  selected  as  one  of  c he  main  tools  in  the  assessment 

procedures  (Ref,  5-7).  The  dac ign  of  an  automated  sequential  filter 
sampler  was  completed  end  a  templing  syst  m  for  sequential  programmer 
collection  wss  initiated,  and  is  discuas»d  in  the  body  of  this  report. 

The  seseesment  procedure  if  based  upon  stirred  sedimentation,  Therefore 
theoreticsl  and  experimental  work  wee  performed  to  determine  the  effac- 
tiveneea  of  chamber  stirring  by  fans  (Ref.  4,5).  Analytic;)]  studies 
wars  performed  ktt  permit  scaling  to  large  test  chambers  of  the  results 
of  mixing  experiments  performed  in  the  small  laboratory  chamber  (Ref.  7). 

Analytical  end  experimental  work  was  Initiated  also  on  the  effects 
of  aleccrical  charging  of  aeroaol  particles  and  their  behavior  in  test 
chambers  -  with  electrically  conductive  and  grounded  walls  (Ref.  7).  This 
study  was  aimed  specifically  at  assessing  the  influence  of  particle 
qharge  on  the  validity  of  the  mass  decay  analysis  for  the  evaluation  of 
explosively  generated  aerosols.  The  initial  phases  of  this  werk  ;erp 
directed  to  establish  an  analytical  model  for  the  combined  de  n‘  on  ef¬ 
fects  produced  by  space  charge  induced  fields  and  gravitc; ionK  settling. 
Experimental  activities  were  concentrated  on  space  charg  measurements 
by  means  of  a  field  mill  and  its  correlation  with  the  mass  decay  of  bath 
intentionally  charged  and  uncharged  aerosols  in  a  small  test  chamber. 

A  review  of  the  present  state  of  knowledge  with  respect  to  the 
deposition  properties  of  aerosols  in  the  respiratory  tract  was  performed 
and  an  extensive  bibliography  on  the  subject  was  assembled,  This  liter¬ 
ature  analysis  pointed  out  various  models  which  have  been  suggested  for 
particulate  deposition  profiles  within  the  human  body  and  the  inherent 
complexity  of  the  biological  structure.  The  pertinence  of  aerodynamic 
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■  v.1?  discussed  in  tt.  ms  of  the  various  deposition  mochai'ism*  .mu  ;  n< 

ti  L  r  t  significance  of  sedimentation  and  impaction  methods  was  r*s  tab- 
(Lcf.  7). 

7 1  ie  ti.'Lu  produced  by  the  .issossmont  system  must  bi  usah  ie  as  input 
to  .i  rviteri.il  balance,  A  series  of  tests  was  performed  in  t  ini  large 
’creed  tret  chamber  ir,  wh  ich.  a  variety  of  tecimi^uos  were  etvimtrd  l  e. 
several  attempts  to  obtain  a  material  balance.  A  workable  method  was 
determined . 

Raoid,  reliable  chemical  analytical  methods  are  the  cornerstone  of 
the  assessment  procedures.  Several  areas  related  to  analysis  have  been 
examined,  Preliminary  work  to  apply  flame  pnotometry  as  an  analytical 
procedure  for  the  determination  of  total  phosphorus  is  reported  herein. 

A  new  method  Co  replace  the  oxidation  seep  in  the  current  phosphorous 
procedure  was  evaluated.  The  entire  filter  and  deposition  material 
including  carbonaceous  material  is  digested  in  a  perchloric  acid  mixture 
to  assure  that  all  the  phosphorus  is  assayed.  This  method  is  also  dis¬ 
cussed  later.  Finally,  the  Technicon  AutoAnalyzer  system  was  integrated 
into  the  chemical  evaluation  chain  and  procedures  were  developed  for  the 
detection  of  various  dissemination  materials  (Ref.  5). 

Laboratory  teating  of  candidate  devices  led  to  work  in  two  areas 
related  to. aerosol  generation.  Preliminary  experimental  work  was  per¬ 
formed  using  saccharin  and  subsequently  DOP  aerosols.  Several  generating 
procedures  and  apperatus  were .evaluated  to  produce  aerosols  in  the  size 
region  of  interest  (a  few  microns) .  The  condensation  generator  was  re¬ 
jected  because  of  an  insufficient  rate  of  generation.  A  cascade  impactor- 
spray  nozzle  combination  was  tested  and  found  to  produce  reliable  aerosols 
but  its  use  was  discontinued  because  an  excessive  fraction  of  particles 
was  smaller  than  1  micron.  Finally  DOP  uerosolization  was  accomplished  ' 
using  pressurized  dissemination  with  Freon-12.  These  tests  indicate  that 
particle  size  properties  may  be  controlled  by  varying  the  relative  amounts 
of  DOP  and  Freon-12  (Ref.  5). 

The  problem  of  unwanted  counts  registered  by  the  Royco  Light  Scat¬ 
tering  Particle  Analyzer  with  latex  calibration  spheres  was  studied. 

The  evidence  points  to  a  large  spurious  population  of  particles  smaller 
chan  0.1  micron  which  apparently  is  capable  of  triggaring  the  circuitry 
of  the  instrument.  Selective  filtration  appears  to  promise  a  solution 
to  this  problem  (Ref.  5) .  „ 

Finally,  a  new  application  of  beta  absorption  techniques  for  mass 
analysis  is  being  investigated  and  promises  to  provide  useful  backup 
information  for  the  mass  measurements  made  by  filtration  techniques. 
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1 J  ! .  DF-SCR I P HON  OF  ACT  IV  i  TICS 


A.  AUTOMATED  SEQUENTIAL  F1I.TF.R  SAMPLING 


d  M  L  Ulil  J  L 


The  design  was,  completed  and  assembly  was  initiated  oi  an 
rillv.  system  for  sequential  programmed  collection 


Man  ml  sampling  procedures,  although  adequate  for  experimental  purposes, 
are  not  to  be  preferred  for  routine  and  repetitive  measurements.  Dif¬ 
ficulties  are  especially  appa^cnc  daring  the  initial  stages  of  sampling; 
e.g.,  shortly  after  explosive  aerosol  generation,  when  a  maximum  of  in¬ 
formation  has  to  be  gathered  in  a  minimum  of  time. 


lhe  sampling  system  has  been  designed  to  operate  as  follows; 

(!'  Preloaded  filter  holders  in  a  special  container  will  be 
introduced  through  an  appropriate  opening  into  a  sealed  "glove  box"  at¬ 
tached  to  the  sice  of  the  assessment  chamber.  This  glove  box  will  con¬ 
tain  a  multiple  filter  sampling  head,  mounted  on  a  base  plate  riding  on 
slide  tracks  The  loaded  filter  holders  will  be  inserted  ,.nto  the 
sampling  head.  A  door  between  the  glove  box  and  the  chamber  will  be 
closed  during  the  dissemination  event  to  protect  the  sampler  from  shock 
waves,  shrapnel,  and  ocher  disturbances  associated  with  dissemination. 


(2)  .  Prior  to  dissemination,  a  manually  operated  switch  will 
start  the  sampling  pump  and  apply  voltage  to  electronic  circuitry  as¬ 
sociated  with  the  sampling  operation  (fan  tachometers  ,  flowmeter  circuitry, 
etc . ) . 


(?)  About  0.5  minutes  after  time  zero,  an  operator  will  open  the 
door  separating  the  glove  box  from  the  chamber  and  slide  ouc  the  sampling 
Cray  with  the  multiple  filter  head.  Thia  action  will  trigger  the  starting 
of  a  sequence  timer  which  controls  the  operation  of  the  solenoid  valves. 
Each  filter  will  operate  in  sequence  and  sampling  will  stop  automatics lly 
at  the  completion  of  the  run  (about  30  minutes  later)  . 

(4)  The  dissemination  event  will  initiate  the  operation  of  a  print 
out  timer-flow  counter  to  correlate  individual  filter  sampling  with  time 
after  aerosol  disarm! nat ton. 

(5)  The  sampling  tray  will  then  1  retracted  manually  into  the 
glove  box ,  the  filter  holders  removed  from  the  sealed  environment  and 
placed  into  a  special  container  and  removed  from  the  bottom  of  the 
glove  box, 

The  sampling  sequence  will  be  programmed  by  means  of  plug-ina  inserted 
into  the  front  panel.  This  feature  allows  the  choice  of  any  number  of 
different  sampLing  sequences  such  that  the  12  filter  samples  with  indi¬ 
vidual  duration  of  1/2  minute  can  be  distributed  in  any  way  over  a  total 
sampling  period  of  up  to  50  minutes.  In  addition  the  individual  sampling 
period  can  be  extended  to  1  minute  and  the  total  to  100  minutes. 
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,\  routing  aver  with  automatic  indexing  will  be  in  front  of 
the  filter  heuus  to  protect  the  non-operating  filters  from  direct  exposure 
to  the  jcrosol  cloud.  During  the  intermediate  dead  time  (when  no  filter 
ia  sampling)  this  protective  disc  will  be  covering  all  filters.  The 

1  <  na  ^  «»•  »  e  e  a»i«  »■  S  **  *«  *•  A  w»  <■»  *4  W«»  ►  W  «  r  .*%»►<  a  <  i  1  .a  *•  r-  m  n  W  mx 

- - 0  •  •-  *  —  *  *  w  r  -  ~a  ~  ■“  ~  j  -  *  • w  r-^  -  w  v****r  -■*■“© 

sequence  chosen. 

Previous  experiments  have  established  quantitatively  the  fol¬ 
lowing  sampling  parameters r 

(1)  Total  sampling  period,  after  time  zero,  of  approximately 
30  minutes . 

(2)  Sampling  intervals  of  0.5  minutes. 

(3)  Minimum  desirable  number  of  sampling  intervals,  ap¬ 
proximately  12. 

(4)  Initiation  of  first  sampling  interval,  about  1  minute 
after  time  zero  tc  allow  for  adequate  chamber  mixing. 

(5)  Initial  sampling  repetition  rate,  one  every  minute, 
decreasing  with  time. 

(6)  Sample  flow  rate,  approximately  1  to  2  cfm  for  a  1-1/2  inch 
diamatsr  filter. 

Several  alternative  valvins  mechanisms  for  the>  automated  se¬ 
quential  sampling  system  have  been  explored.  The  Gelman  sequential 
sampling  head  (type  No.  1)  consists  Of  a  circular  array  of  12  needle 
values  actuated  by  a  "geneva1’  mechanism.  This  device  was  rejected,  among 
other  reasons,  because  it  presented  an  excessive  resistance  to  flow, 
making  it  impoaaible  to  obtain  the  desired  flow  rate.  A  flow  rate  of 
only  0.74  cfm  required  a  preisure  irop  of  23  inches  of  acroaa  the 
sequential  valve.  No  other  adequate  sequential  valve  mechanism  seemed 
to  be  evailable,  and  thus  s  system  of  timer  actuated  solenoid  valves  was 
chosen.  Flow  measurements  were  performed  on  severs!  solenoid  valves  to 
establish  their  f.'ow-pressure  drop  characteristics,  end  a  valve  with  a 
5/16-inth  orifice  vith  a  pressure  drop  of  3.6  inches  of  Hg  st  2  cfm  was 
chosen . 


It  w  •%  decided  to  measure  the  total  air  volume  sampled  througft 
each  filter  or.  _n  automated  !.  *ei«  by  ...nan  a  of  a  turbine  flowmeter  and 
to  provide  electromech  .nical  counter-recorder  display  of  the  total  flow. 
Sample  concentration  vi’i  be  normv1,zed  on  the  basis  of  the  in1  grated 
flow  through  each  filter.  Au^omate^  constant  flow  rate  by  feedback  con¬ 
trol  appears  unnecessary  and  would  add  unwarranted  complexity  to  the 
ays  tern. 
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Although  only  one  system  of  this  type  is  being  contemplated 
ut  this  point,  any  number  of  them  could  bs  included  for  simultaneous 
sampling  at  various  sices  within  the  chamber.  The  entire  sampling 
system  describee  above  is  designed  for  maximum  reliability,  sturdiness 
and  simplicity  of  operation,  making  maximum  use  cf  commercially  avail¬ 
able  devices  ,  wherever  possible, 

o  \t  k  Trp**t  nv  wv-n^v*  r\ir:-%rmr> 

L>  i  riAlwhUlia  uAMftkiCLi  i-  A-i  1.U  ii'ICkUJ 

A  series  of  tests  wss  performed  in  the  large  aerosol  test  cham¬ 
ber  in  which  an  attempt  was  made  to  recover  or  eccount  for  all  of  the 
bomb  fill  material  by  sampling  tha  chamber  floor,  walls  and  fans  as  well 
as  the  aerosol. 

The  fill  materiel  used  for  most  of  these  tests  was  a  1  percent 
water  solution  of  fluorescent  uranine  dye.  At  first  the  floor  was 
sampled  by  covering  it  with  trays  which  were  separately  washed  for  re¬ 
covery  of  the  dye.  After  the  aerosol  was  sampled,  the  sum  of  the  floor 
and  aerosol  recoveries  was  less  than  100  percent  and  was  variable  from 
test  to  test.  Table  1  summarizes  the  results  of  these  tests. 

Another  procedure  was  tried  which  avoided  the  tedious  work  of 
washing  about  300  traye. 

(1)  A  bomb  was  fired  end  the  aerosol  sampled. 

(2)  The  chamber  floor  was  flooded  with  several  hundred  gallons 
of  water.  After  thorough  mixing  to  dissolve  tha  settled  dye  ths  solution 
wet  eamplad. 

(3)  A  spray  was  used  to  waah  the  fens  with  an  additional  meaa- 
urtd  amount  of  water  which  was  mixed  into  ths  watar  already  on  tha  floor. 
Another  sample  was  taken. 

(4)  The  overhead  rotary  noexle  was  turnad  on  to  wash  tha 
celling  and  walla  of  the  chamber  into  the  previous  washings  and  a  third 
sample  was  taken. 

(5)  A  carefully  measured  amount  of  the  bomb  fill  solution  was 
added  to  the  washings  which  were  mixed  and  sampled  again  to  provide  an 
internal  standard  for  the  rest  of  the  analyasa. 

(6)  Tha  aerosol  was  again  sampled. 
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FLOOR  SAMPLED  WITH  TRAYS 


Up  to  this  poinc  the  chamber  was  not  vented,  and  all  the  dye 
In  the  boob  plus  the  known  added  amount  was  assumed  to  be  in  the  water 
accumulated  on  the  floor  except  for  a  very  small  amount  still  suspended 
in  the  air. 

Table  2  gives  the  results  of  three  tests  using  this  sequential 
washing  technique.  It  appeared  that  for  the  high  mass  ratio  31  type  boob 
substantially  sll  the  dye  was  recovered.  With  the  low  mass  ratio  No.  23 
sphere,  less  than  40  percent  of  the  fill  was  accounted  for.  It  was  con* 
eluded  that  tha  dy«  was  decomposed  or  dtmagsd  by  the  dissemination  in 
this  bomb. 

A  20  percent  solution  of  phosphoric  acid  in  water  was  dis¬ 
seminated  in  testa  457  and  459.  Sample*  were  analyzed  on  the  Technicon 
for  total  phoaphorua.  Again  the  total  recovery  was  low.  In  this  esse 
we  suspected  that  the  acid  raactad  with  the  steel  floor  of  the  chamber 
to  form  Insoluble  compounds  which  adhered  to  the  floor. 

Prior  to  these  tests  and  on  another  program  a  few  tests  were 
run  with  Bis  and  VX  filled  bembs  wherein  after  the  aerosol  was  sampled 

and  wanted  from  the  chamber  only  the  floor  was  waahed  with  about  10  gal- 
lens  of  isopropanol.  The  washing  solution  was  sampled  and  analyzed  In 
the  same  manner  as  tha  filters.  Total  r sc ovaries  were  eigniflcently 
lower  than  100  percent,  suggesting  a  significant  wall  loss.  Tha  later 
teats  with  uranine  solution  have  indicated  that  a  significant  and  vari¬ 
able  fraction  of  the  fill  materiel  anda  up  on  the  walls  and  fane. 

Tha  sequential  washing  of  the  floor,  fans,  and  walla  followed 
by  sampling  of  the  air  Is  necessary  for  a  good  material  balanca.  This 
procedure  is  workable  but  presents  the  following  problems. 

(1)  A  dense  end  lethal  serosol  Is  encountered  in  carrying  out 
this  procedure  manually.  To  avoid  this »  seperete,  remotely  controlled 
weshlng  fecllltles  must  be  provided  for  floor  and  fans. 

(2)  Large  amounts  of  solvent  ere  needed.  A  minimum  of  100 
gallons  par  test  ia  estimated  for  the  large  chamber.  It  is  desirable 
that  thia  aolvent  ba  cheap  and  not  fora  an  explosive  mixture  with  the 
chamber  sir,  end  suet  b«  comps t Ibis  with  the  chemical  analysis.  Tha 
conventional  analytical  solvent,  isopropenol ,  ia  unsatisfactory  because 
of  expanse  and  for  reasons  of  safety. 

Water  and  aqueous  solutions  of  hexylene  glycol  and  bleach 
are  being  considered  as  inexpensive  and  seft  substitutes  for  the  alcohol. 
Experiments  were  conducted  to  evaluate  the  solubility  of  VX  in  aolutiona 
of  various  strengths  and  the  ability  of  the  analytical  method*  to  recover 
VX  in  each  of  the  solutions.  The  analytical  methods  of  particular  interest 
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were  the  Si  hoc- re mann  me L hod  for  active  agent  and  the  ferric  perchlorate 
hydrogen  peroxide  method  for  total  phosphorous.  Results  are  shown  in 
Figure  1.  VX  Is  stable  in  a  solution  of  50  percent  glycol  even  ii  stored 
over  u  weekend,  and  water  alone  is  an  adequate  solvent  if  analyses  can  be 
perfe  mod  itnmed  is  Ce  ly  . 

VX  was  dissolved  in  a  hypochlorite  decontaminating  solution 
and  was  analyzed  by  the  total  phosphorous  method.  One  hundred  percent 
of  r ha  vv  w.ii  recovered,  indicating  that  the  "drccr.11  conld  be  need  if 
only  phosphorus  balance  was  required. 

C.  PERCHLORIC  ACID  OXIDATION  0?  FILTERS  FOR  TOTAL  PHOSPHORUS 

One  possible  explanation  for  the  difficulty  in  achieving  a 
material  balance  in  chamber  tests  is  poor  sxtraction  efficiencies  because 
of  tenacious  retention  of  agent  and  agent  decomposition  products  on.  the 
carbonaceous  decay  products  collected  on  filters  during  sampling.  Since 
in  the  normal  analytical  procedure  only  the  extract  it  analysed,  the 
retained  material  would  not  be  accounted  for.  A  new  method  of  oxidation 
was  suggested  by  Mr.  A.  Koblin  of  Edgewood  Arsenal  in  which  the  entire 
filter  ia  digested  in  a  perchloric  acid  mixture.  The  carbon  i*  oxidized, 
and  thus  releases  the  suspected  absorbed  materials  .  The  laboratory 
procedure  for  the  method  wee  established  for  bis  end  then  et  the  suggestion 
of  Hr.  Koblin  for  trioctylphoephat*  (TOP).  Firat  determinations  were 
made  on  solutions  of  Bit  and  later  with  TOP  with  excellent  results  in 
both  cates  ns  shown  in  Figure  2.  The  materials  war*  then  applied  to  clean 
filters  and  carbon  coated  filtara  and  the  ability  to  recover  the  material 
was  evaluated.  Again,  results  wars  excellent  with  average  Bia  recoveries 
of  W  portent  in  12  samples  with  an  res  deviation  of  1  percent  end  an 
average  TOP  recovery  of  97  percent  from  7  samples  with  an  rate  deviation 
of  1  percent. 

Following  the  eucceae  of  this  method  for  Bis  and  for  trioctyl* 
phosphata  (TOP)  its  application  was  attempted  on  VX,  The  following  kinds 
of  temples  were  takent 

(1)  VX,  weighed  into  ■  teat  tube, 

(2)  VX,  weighed  onto  a  filter  paper,  and 

(3)  VX,  weighed  oneo  a  carbon  coated  filter  paper, 

In  all  cases,  the  temple  wet  covered  with  the  digestion  mixture 
and  oxidized.  Six  aliquots  of  each  ware  tiken,  ranging  in  site  from  25a 
to  1.00  ml  and  the  color  reaction  completed.  The  results  ware  calculated 
from  the  best  values  of  the  absorbance,  that  is,  when  the  value  was  be¬ 
tween  0.050  and  0.900.  Some  of  the  results  are  shown  in  Table  3.  There 
is  much  to  wide  a  range  of  recoveries  for  the  method  to  be  useful  to  us, 
unless  the  performance  can  be  improved. 
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Figure  1.  Stability  of  IX  la  aqueous  solutions  of  hexylene  glycol. 
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Since  L he  method  works  so  well  for  Blu  and  for  TOP  the  dif¬ 
ficulty  has  Co  be  In  the  oxidation  step.  Communication  with  Adgewoou 
Arsenal  revealed  similar  difficulties  in  applying  the  perchloric  acid 
oxidation  to  VX.  A  literature  search  of  the  procedure  was  made  which 
included  3  search  of  the  history  of  the  phosphomoiybdace  color  reaction 
and  the  materials  which  interfer.  Several  possible  interferences  were 
uncovered,  however,  it  was  concluded  thac  the  problem  is  not  one  which 
can  be  solved  simply.  The  analytical  chemistry  aspects  of  Luis  piugrtUu 
permit  application  of  methods  but  do  not  allow  for  the  development  of 
new  analytical  methods-  It  is  suggested  that  the  non-ruproducibil  ity 
of  the  pcr'ehioric  acid  oxidation  of  VX  be  investigated  in  the  Edgovood 
Arsenal  laboratories.  It  is  &  desirable  technique  for  assessment  which 
we  would  be  pleased  to  apply  if  it  is  worked  out  ir.  rime. 

D.  AEROSOL  MASS  EVALUATION)  BY  BETA  ABSORFTION 

Experiments  were  reported  in  the  First  Annual  Report  (Ret.  4) 
on  the  application  of  beta  absorption  measurements  for  the  evaluation  of 
aerosol  mass  collection  on  filter  substrates.  These  tests  appeared  to 
be  discouraging  in  term3  of  the  sensitivity  required  for  typical  sampling 
volumes . 

The  technique  was  subsequently  reevaluated  in  the  direction  of 
u*4njt  a  low  sbeorption  collection  substrate  end  high  collection  concen¬ 
tration  on  a  small  area. 

Initial  experimental  date  were  taken  to  determine  the  feasibility 
of  using  uete  radiatirn  absorption  through  an  eeroeol  deposit  collected 
by  impaction  n  a  v«  /  thin  mylar  substrate,  The  thin  mylar  bate  (1/4 
mil  '•hick)  w<i  ueed  to  minimite  the  ebaorption  due  to  the  substrate, end 
improve  the  seneitivity  of  thia  method.  Using  a  1  inch  diameter  C1*  beta 
source  end  collecting  treated  saccharin  aeroeol  on  mylar  diacs  mounted 
on  the  impaction  plates  of  a  GCA  cascade  impector,  it  was  found  thet  thia 
method  showed  an  improvement  over  the  previously  reported  experiment. 

The  improvement  in  aenuitivity  etill  provided  only  marginal  feasibility, 
but  che  results  pointed  towards  modif ieetlone  in  source-detector  geometry 
which  would  provide  the  desired  sensitivity  for  the  measurement  of  mass 
of  deposited  aeroeol.  The  eiaentlel  modification  appeared  to  be  a  signif¬ 
icant  decrease  of  the  area  of  the  source  in  order  to  decrease  the  amount 
of  particulate  deposition  to  e  practical  level.  A  2  microcurie  C  source 
with  an  active  diameter  of  lass  than  1/4  inch  was  acquired. 

It  was  found  necessary  to  include  eome  collimation  of  the 
r  idiation  into  the  Geiger  detector  tube  to  prevent  errors  due  to  electron 
scattering.  The  collimator  consisted  of  a  short  cylinder  (0-8  cm  long) 
with  an  inner  diameter  of  about  0.55  cm. 
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The  calibration  curve  obtained  with  the  above  arrangement, 
us  in,,  mylar  films,  is  nhovn  on  the  accompanying  graph  (Figure  1).  The 
equation  of  this  curve  is  TR  ■  100  exp  (-0.2156),  where  Tr  is  the  rel¬ 
ative  tranr.miasion  in  percent,  and  5  is  the  ares  density  of  the  absorber 
in  mg/cm-.  Gnat Cenuated  counting  rates  were  about  A3  ,000  counts  per 
minute  which  for  a  1  minute  period  corresponds  to  a  total  error  of 
2\A3,GGu  *  414  wiiiwli  corresponds  to  a  relative  transm***."!  nn  of  ab.-ut 
99  percent,  thus  s  2  percent  reduction  of-  transmission  appears  to  be 
an  approximate  practice!  limit  of  detection.  Transmission  of  98  percent 
corresponds  to  an  area  density  of  about  0.1  iag/cm^  which  in  turn  Cor¬ 
responds  to  20  micrograms  of  materials  deposited  on  a  0.2  car  area  (active 
source  area),  This  sensitivity  would  be  somewhat  degraded  by  the  presence 
of  the  mylar  substrate  or  collection  surface.  For  a  6,3  micron  (1/4  mil) 
mylar  film  it  is  estimated  that  the  limit  of  detectability  would  be  be¬ 
tween  30  and  40  micrograioe  on  a  0.5  cm  diameter  area.  This  result  scemc 
quite  encouraging.  The  next  step  should  be  to  collect  aerosol  or.  mylar 
film  and  corroborate  the  above  conclusions. 

The  seusicivicy  obtained  w.!  th  the  above  geometry  corresponds 
to  ;•  minimum  of  about  1.0  to  1.5  n.icrograms/litcr  (flow  rate  30  1/m  for 

1  minute  sampling  time)  compered  to  the  figure  of  167  miccograma  por 
liter  cited  earlier  (Rc-f.  4).  S*ncc  the  estimated  minimum  desirable  l 
sensitivity  .was  about  5  microgroma / lit-,  r ,  this  technique  a  asms  to  offer 
great  promise  and  will 'be  pursued. 

K.  F.I.EC7RCSTAT1CS 

1.  Experimental 

Tito  experiments  initiated  during  the  seventh  quarter  on 
space  charge  decay  measurements  in  the  small  test  chamber  were  continued 
with  the  main  emphasis  on  correlating  those  determinations  with  “urosol 
mass  decoy  measurements .  The  prime,  objective  of  these  tes  ts  was  to 
establish  the  influence  of  electric  particle  charges  on  the  depletion 
rate  of  a  confined  aerosol  cloud  and  thus  to  ddtermine  the  validity  of 
the  mass  decay  method  for  the  assessment  of  aerosols.  Experience  at 
lidgewood  Arsenal  and  GCA  Corporation  suggests  that  anoraaloue  behavior 
is  not  encountered  when  working  with  explosively  oisseminatad  liquidB 
but  is  ouch  more  likely  for  the  case  of  pneumatically  diseerainated  pro- 
sired  solids.  Consequently,  tha'  latter  was  selected  for  study, 

The  first  set  of  experiments  consisted  of  the  pulsed  in¬ 
jection  of  0.25  gram  samples  of  saccharin  into  a  small  cylindrical  test 
chamber  (50  cm  diameter  x  60-8  cm  height,  of  galvanised  steel).  The 
aerosol  was  injocted  through  a  corona  charging  section  (Figure  4). 
Concurre.it  measurement  of  apace  charge  with  the  field  mill,  and  filter 
samples  for  mass  decay,  were  taken  over  a  period  of  up  to  10  minutes 
after  the  aerosol  injection  time. 


Figure  4.  Needle  point  charger. 


The  results  are  shown  in  Figures  5  and  6.  The  net  space 
charge  measurement  was  made  with  a  field  mill  located  on  the  axis  and 
at  the  top  of  the  chamber.  For  this  geometry,  the  net  space  charge 
density  (c^)  is  given  by 

2E  e 

Co  *  ~ P  *  (nk,} 


Where  £  is  the  field  strength  measured  by  the  field  mill,  l  is  the 
length  of  the  chamber  and  e  la  the  permittivity  of  free  space.  The 
maa^  concentrations  were  determined  by  filter  samplee  (2.5  liter  volume) 
collected  at  1  minute  interval*.  In  figures  5  and  6,  the  circled  numbers 
ere  the  run  nuabors  and  the  numbers  beside  the  circles  are  the  voltage* 
applied  to  the  high  voltage  electrode  of  the  charging  section.  The  sign 
of  po  is  the  same  at  the  sign  of  the  high  voltage  electrode. 


It  can  be  seen  in  both  Figurea  5  and  6  that  the  data  are 
split  into  two  groups.  Since  there  seemed  to  be  no  obvious  explanation 
for  this  effect,  the  experiment  was  rerun  to  determine  if  the  effect  was 
real  or  coinc idet.tal .  The  results  of  thi*  experiment  are  aho-r,  in  Figure 
7.  In  this  case  no  significant  correlation  between  mass  decay  character- 
...4a.  tic  a  aari  charging  wa»  ebaerved.  Since  the  above  set  of  results  were 


inconsistent  end  inconclusive  a  new  improved  experiment  was  performed,. 
The  new  sxperiaant  cowiitW  ef  .  the  direct  meaaureaeot  »1  tha  aafoeol 
lost  to  the  cheats r  floor  and  walla  and  tha  comparison  of  the  results 
for  casta  where  there  is  deliberate  charging  of  the  aerosol  with  caese 
where  there  it  no  delibtretc  charging.  If  gravity  is  the  main  cbtvtzl-  ; 
butter  to  ehe  loa*  process »  the  bulk  of  the  lost  aerosol  should  be  found 
on  tlie  “ftfcfdir *’•••;. tt  eTeetrdgti tied  ere  important  is  the  loss  process ,  this 
should  be  reflected  in  high  wall  losses. 


In  order  to  measure  the  floor  and  wall  deposition  of  sac* 
chatin,  1/32  inch  thick  aluminum  plates  were  placed  on  tjhe  floor-  end  at¬ 
tached  to  tha  wall  of  the  chamber.  The  operating  geometry  of  the  plates 
in  the  chaster  is  depicted  in  Figure  8. 


It  wee  suggested  thee  the  charging  aection  used  (Figure  4) 
might  ba  ambiguous  as  to  where  the  corona  discharge  takes  place.  In 
order  to  avoid  this  difficulty  a  cor^entric  charger  as  shown  lu  Figure  9 
waa  constructed.  In  addition  it  was  found  preferable  to  monitor  the 
charging  current  instead  of  the  voltage. 
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Figure  8.  Collection  plate  geometry. 


The  procedure  for  these  deposition  tests  wss  as  follows. 

About  0,25  grams  of  saccharin  were  injected  into  the  chamber  through  the 
enneentrie  charging  section.  For  five  runs  the  charger  was  not  orcr.itcd, 
while  for  five  ocher  runs  the  center  electrode  was  positively  bi.ispd 
•uch  chat  a  corona  current  of  10  microamps  was  measured  flowing  to  ground. 
After  injection,  the  charging  section  was  turned  off  (for  the  five  runs 
where  applicable)  and  the  saccharin  was  allowed  to  settle  for  10  minutes 
with  the  mixing  fan  running.  At  thia  point,  the  deposition  places  were 
removed  and  the  mats  on  the  wall  plates  and  the  mass  on  the  floor  plates 
were  measured.  Tht  results  of  these  experiments  are  shown  in  Table  4. 

Neglecting  the  eighth  run,  which  shows  a  somewhat  anom¬ 
alous  behavior,  the  following  conclusions  are  indicated.  Wall  accumu¬ 
lations  are  roughly  twice  as  great  when  the  particles  are  artifically 
(rrrons)  charged,  indicating  pronounced  electroatatlc  deposition.  Floor 
mass  densities  are  only  slightly  affected  by  artificial  charging. 

The  effect  of  particle  charge  on  a  volume  sampling  assess¬ 
ment  ayjtem  is  given  by  total  mass  loss  rather  than  mass  density  depo¬ 
sition.  Thus,  averaging  the  first  five  end  cite  last  four  (omitting  the 
eighth  run)  runs  And  multiplying  by  the  appropriate  floor,  calling. and 
wall  areas,  one  obtains  ths  reaults  ahown  in  Table  5.  The  objective  of 
the  experiment  was  to  look  at  the  wall  and  floor  depositions  and  to  try 
to' observe. a  difference  between  the  deliberately  charged  and  the  naturally 
charged  cages.  Unfortunately  the  cloud  concentration,  tha  amount  of 
material  on  tha  fane,  and  tha  amount  of  residue  in  the  dlaaemlnator  and 
charging  section  were  not  evaluated.  Thia  was  ragratsble  aincs  the  a- 
mounee  ahown  in  Table  5  do  not  add  up  to  tha  ISO  og  original  aample  weight'. 
Teats  should  be  conducted  in  which  all  the  particles  are  completely 
accounted  for. 

Nontthalsaa,  a  vary  aignlflcant  difference  in  deposition 
pattern  via  observed  for  tljs  cast  of  deliberate  charging.  It  la  claar 
that  for  the  case  of  pneumatically  disseminated  preaired  povdare  the 
electrostatic  forces  art  pi#  ing  a  definite  role  in  removing  particles 
from  the  test  volume.  Further,  theae  raaulta  are  averaged  over  the  time 
necessary  to  remove  oor.t  of  tha  particles  from  the  volume.  At  early  times, 
tha  effects  of  electrostatic  forces  would  be  even  more  pronounced. 

No  electrostatic  meaeursoenta  have  yat  been  made  in  the 
large  test  chamber.  Sinca  the  dimensionlese  ratio  r^  qZL/e  og  is,  for 
constant  charge  density,  five  '•imas  as  large  for  the  large  fank  as  for 
the  email  one,  the  indications  are  that  electrostatic  forces  will  play 
some  role  and  may  even  dominate  thf-  assessment  in  the  large  tank.  The 
significance  of  this  parameter  is  explained  below  in  the  anelytic  treatment. 


TABLE  5 


MASS  DEPOSITED  ON  FLOOR  AND  WALLS  FOR  TWO  CHARGING  CONDITIONS. 


Total  Mass  Deposited 

Total  Mass  Deposited 

on  Wall 

on  Floor 

Intentionally  (corona) 
Charged  Particle* 

57  mg 

63  mg 

"Naturally"  (no  corona) 
Charged  Particle* 

26  mg 

63  mg 

Tha  above  considerations  do  include  a  number  of  tenuoua 
aaatnnptiona .  Thu« ,  it  would  be  neceatary  to  perform  electrostatic 
measurement*  in  tha  large  chamber  to  teat  theae  conalderatlona  directly. 

A  conclusion  drawn  from  theae  teata  la  that  for  the  caae  of  pneumatically 
diaaaalnated  dry  powder#  a  aerioua  error  could  be  introduced  if  only 
gravitational  eettllng  la  aaauaed  for  aaaa  concentration  decay  without 
a  correction  procedure  for  loaa  due  to  electrostatic  effects  on  the  wall, 
calling,  and  floor  of  tho  chamber. 
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Future  work  in  this  area  is  expected  to  clarify  thic 
problorj,  and  tests  along  these  linos  are  being  planned  for  the  large 
(213  n  )  chamber . 


2 .  Analy ticsl 

A  theoretical  treat-rent  of  the  settling  of  particles  by 
combined  gravitational-electrostatic  has  been  made. 

a.  Monooolsr  Charging 

The  probl-m  being  considered  here  is  the  following. 


For  the  case  where  an  oeroanl  cloud  of  particulate 
matter  is  introduced  into  a  container  having  an  electrically  conducting 
inner  surface,  it  is  desired  to  obtain  th-*  law  by  which  particles  arc 
removed  from  the  volume  and  deposited  on  the  aurfuce.  It  will  be  assumed 
that  the  interior  of  the  container  is  divided  into  two  regions: 


(11  A  thin  boundary  layer  In  which  the  gee  ir,  the 
container,  is  assumed  to  be  cwtionleie. 


(2)  A  region  outside  the  boundary  layer  In  vhlch 
the  gas  is  stirred  in  such  a  way  sc  to  create  a  homogeneous  distribution 
of  particulate  matter .  ' 

It  may  be  remarked  at  chit  point  that  this  division 
may,  as  far  as  practice  is  concerned.  ba  open  to  some  question.  Beside 
questions  concerning  lha  rsalii ability  of  homogeneity  In  the  Interior 
region,  there  1*  also  «  problem  concerning  the  relative  thicknese  of 
the  boundary  layer  end  a  "pars is ttnre  length"  of  a  particls  once  set  into 
motion.  Ih  order  that  thia  division  be  valid,  the  boundary  layer  must 
be  many  tiraee  the  diatenoe  over  which  particle*  arc  brought  to  a  conatant 
velocity.  For  present  purpoiaa,  thee*  questions  will  not  b«  examined 
further,  and  thie  division  will  be  mad*. 

Coneldar  fir at  the  arrangement  ehown  in  Figure  10, 
a  cylindrical  container  end  particle*  of  a  single  else  and  charge.  The 
more  complicated  ceiae  of  particle*  charged  both  positively  and  negatively 
and  of  varying  magnitude  of  charge  and  diameter  will  be  deferred  to  a 
later  report.  The  ioll<ntng  rotation  supplement*  the  notation  glvth  in 
Figure  10: 

R  -  Cylinder  radlua 

!.  -  Cylinder  height 

A.  -  Cylinder  ...d  area 

Aj  -  Cylinder  aide  „r#* 

N*  -  Total  inutantentoua  number  of  undepoalted,  charged 
particle*  in  container 

n  -  ini tanteneoue  number  density  of  undeposited  chsrged 
part  idea  ir  container 
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initial  number  density  of  charged  particles 
charge  per  particle  (assumed  constant) 
particle  diameter 
Particle  mass 

gravitational  acceleration 

coefficient  of  viscosity 

drift  velocity  due  to  gravity 

drift  velocity  due  to  electric  field 

average  value  of  electric  field  at  the  container  wall 

permittivity  of  free  space  (10'^/36n) 


Olitl7-IOI 


2  7TRL 


Figure  10.  Definition  of  symbola. 


Since  che  p-irticles  are  considered  to  bs  charged,  in 
addition  to  a  downward  gravitational  drift  over  the  lower  area  A^ ,  the 
particles  are  acted  on  by  an  electric  field,  the  source  of  which  is  the 
not  charge  on  the  particles  themselves. 

The  balance  equation  for  deposition  will  be 


d  N  ,  _  d  n  .  < , .  .  .  . 

Tz  m  AiL  Ttm  '  V  vg  *  (2Ai  +  V  n  vf 


CD 


The  drift  velocities  are  obtained  by  assuming  that 
Stoke 's  lav  holds  in  th«  boundary  leytr, 


6n  n  a  v  “mg 
8 


6*  n  a  v. 


qE 


<25 

(3) 


E,  the  value  of  the  electric  field  at  the  surface  of 
the  container  will  for  a  cylinder,  vary  ever  the  surface.  For  our  pur- 
potas  ,  however,  a  repreientative  value'' of  E  mey  be^  obtained  by  eh  ele- 
rjentary  application  of  Game’  law. 


qn  n  R2L  -  (2*  R2  -r  2rt  RL)  E  e 


<*> 


Putting  Iquationa  < ^ ,  (3)  and  (43  into  Bquaci.ua  (4), 
and  replacing  ,  Aj  by  their  equivalanta  in  R,  L,  ons  obtain*  i 


&  .  - 

dt  \in  n  a l. 


'J  “  nT  a. 


■<$>. 


conditions 


The  aolutlon  of  Equation  (5)  aubject  to  the  initial 


is 


n  -n. 


t  a  0 


q2i 

mg  f 


+  1  j  e  ip 


mg 


,pr  n  aL 


n!  q  1 

t  mg 

o 


(6) 
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For  q  **  C  one  obtains 


(7) 


w 


nl q 

1  +  6n  n  a  €  C 

o 


the  ramon  exprmib.'  aaeaciett’d  with  aquere  law  t^saobintcioa.  The 
relative  influence  of  electric  field  and  gravity  e*  early  timer  it  that 
obtJitttd^frcsn  the  dlnenflon&eia  ratio  n.q*L/ag<  u  can  be  aaen  directly 
-by -inepection  ofEquetlc<i  (5>. 

'  iquitien  (6)  applier  only  to  theea  particle*  vUleh  . 

contribute  to  t;  '  net  charge  in  the  tank.  Uncharged  perticlae  ere  removed 
onijr  by  girevitetlonel  settling  and  hence  elveyi  settle  according  to  I 
•lew  of  the  type. given  by  Equation  (7).  tf  the  iniclel  particle  density 
for  neutral  par  tic  lei  la  n5,  than  the  total  number  demitlei  at  time  c 
wilt  be' given  by  * 


b.  Numerical  Example 


Having  derived  equation!  for  the  settling  rate  of 
particle*  under  combined  gravitational-electrostatic  forces,  these 
equations  will  now  be  applied  to  the  settling  occuring  in  two  of  the 
experimental  systems  used  at  GCA,  a  small  cylindrical  metal  tank  and  a 
much  larger  cylindrical  tank,  also  with  an  alectrlcally  conducting 
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interior.  The  examples  trtacsd  are  for  an  aeroiol  consisting  of  sac¬ 
charin  particles,  assumed  to  have  a  single  size  and  mass.  The  pertinent 
numbers  are: 


Small  tank 
Large  tank 


diameter  50  ems ;  height  60  caa 
diameter  6.1  meters;  height  7.3  meters 
5  microns 
1  g/co^ 


Particle  density 
Total  mast  injectad  into  smell  -tank  Q.25g 
Total  meat  injectad  into  large  tank  200g 

3 

Acceleration  of  gravity  ■  10  ega  ■  10  taka 
Single  (spherical)  particle  mess  t  6.25  x  10 


■11 


Coefficient  of  viscosity  of  air  t  1.8  x  10  cgs 


(1)  Exponential  Time  Constants 

T .  i\,n  «.  ?  a  ; 

.  v  v.'.;  *.25  *  Iff'1*  II  103 

■  27  minutes  (small  tank) 

■  6  hours  (large  tank) 

(2)  Dimensionless  Ratio  n^q 

This  ratio  involves  electrostatic  quantities 
und  mk«  units  will  be  used  when  convenient.  '• 

> 

*  3.42  x  lQl®/m*  for  the  small  tank 

■  1,5  x  10^/rn^  for  the  large  tank 

If  it  is  assumed  that  -  1/5,  then 

9  3 

n^  >6.84  x  10  /o  for  the  email  tank 
9  3 

-  3.0  x  10  /a  for  the  large  tank 
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The  average  charge  per  particle,  q,  may  be 
found  approximately  by  a  crude  application  of  Gauss'  lav.  Thus,  if  it 
is  assumed  that  the  field  ia  constant  over  the  whole  interior  surface  of 
the  tank,  Causa'  lav  gives 

i  fR-f-J.)  a  t 

q  .  q - — -2- 


Typicel  valuta  for,E,  found  by  uaa  of  a  field  mill  in  the  email  tank,  era 
of  the  oroar  of  10*  volts/aater,  when  injected  particlta  art  charged  by 
corona  and  of  ehe  order  of  1/10  thia  value  vnan  naturally  charted.  Ut lag 
the  higher  (corona  charged)  value,  ona  obtains  q  ■  1,45  x  10"i°  coulomb 
for  the  email  tank  and  the  large  tank  (tinea  the  charge  dependa  on  the 
charging  procaaa  and  not  the  tank).  Using  all  of  thaae  valuta  one  obtaina 

nr  q2L 

— - — — *  *  15.6  for  the  email  tank 


■  83  for  the  lar^fc  tank 

thus ,  to t  thaee  conditions  a  fraction  n./tu+n, 
of  the  total  injected  particlta  ia  Initially  more  influenced  by  tltctra- 
itsti<  than  by  gravitational,  forces.  Ultimately  of  course ,  with  passage 
of  iim,  Kquation  (6)  revert*  to  an  exponential  law  one*.  aqaln,  but  for 
mudh  tttjrt  r./n*  «  1  a.-#  aaaeaatannt  amta&mn&ca  may  he  difficult  to 
■■■iihe  . 

•  .  •■y '  ■  ■  •  '  •  •,  ’  •'  ■  •• 

c.  lipolar  Charging 

In  thia  taction  the  problem  ^  8Xton&4  to  include 
charges  of  both  aigna.  The  time  nlatoriee  of  tna  charged  particles  are 
doacribad  by  tvo  coupled,  nonlinear,  simultaneous  differentia  V  equations 
for  tb/v  particles  having  plw  and  minus  charge.  £  single  ntegral  ia 
oLtalntd,  and  this  Uada  to  a  decoupling  of  the  aquations;  the  resulting 
equation*  confining  only  a  tingla  dependant  variable  era  still  non¬ 
linear,  and  cannot  be  solved  <si  terms  of  elementary  function*.  Further 
effort  la  required  to  either  solve  these  aquations  numerically,  find  ap¬ 
proximate  methods  of  solution,  or  solve  alaplos  equation*  which  may  ba 
uaad  to  place  <ippe'.:  and  lower  bounds  on  settling  rates.  The  following 
parameters  are  added  to  the  prior  notation  list. 

n  -  number  of  density  of  neutral  particles 

n°  -  number  of  density  cf  positively  charged  particles 

n_  -  number  of  density  of  negatively  charged  particles 

For  further  details  of  the  problem  being  treated,  refer  to  the  treatment 
of  monopolar  cht.g.  Ag. 
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particle*  ar» 


The  equation*  for  deposition  for  the  three  type*  of 


dn 

A1L 

0 

dt 

*  -  A,  n  v 

1  o  g 

(10) 

dn 

A1L 

+ 

dt 

"  *  A1  n+  Vg  ■ 

(2Aj_  + 

V  n+  vg 

(ID 

dn 

r  ° 

if 

v,  >  v  'i 

A1L 

“ 

’  l*V-  <vg  ■ 

•  vf)  if 

h<v|} 

(12) 

For  the  case  v,  ^  v  ,  the  solution  of  Equation  (12)  1*  n  "  constant, 
the  solution  for  n  8reverts  back  to  ona  similar  to  that  ?ound  for  tlia 
case  of  settling  in  the  presence  of  particles  of  ona  sign, 

2 

In  general,  setting  A^  ■  «  R  ,  A„  ■  2«  RL,  using 

Stokes's  law, 


6n  a  •  mg 
6«  a  Vj  ■  qE 

and  Gauaa  law 

q(n+  -  nj  k  %h  •  (In  RJ  +  in  RL)  S  tQ 

ono  obtains 
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,<V‘} 
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and 


no(°)  -  nL 


n  (O')  ■  n . 

f  i 


n  (0)  ■  n. 


Then 


X 

V 

z 

X 


+  rtj  +  n^ 


V"l 


V”i 


n  In. 


n2/ni 

n2/'ni 


•  n3/nt 


lt  q2!. 


n,  q‘L  I 


2  «  '  a*  (R+l)  *  P 
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-X 
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a  .  (° 

dT  L  -z  + 


3  (Y-Z)  >  1 
Z  +  3Z  (Y-Z)  3  (Y-Z)  <  1- 


(10b) 

(Hb) 

(12b) 


The  equation*  are  now  in  fully  non-dimena ionalized  form,  and  the  initial 

condition*  a  re  X«X,Y»Y|Z,2,t"0, 
ooo 


40 


Multiply  Equation  (lib)  by  f3/Y,  Equation  (12b)  by  a/Z  and  add 


e  dY  2  dZ 
Y  d*  +  Z  di 


(13) 


J.11  i.t*£  l  a  l.  mg  diiu  ktitxi  ua  uglug  , 


Y*  Z"  -  C  exp  C-2t)  (14) 
s  a 

Using  the  Initial  conditions,  C  ■  Y  ZQ,  so  that  Equation  (14)  supplies 
a  relationship  between  Y  and  Z.  Usingthis  relationship  uncouples  the 
Y  and  Z  equations. 


dY 

d: 


-Y 


OYk  +  02 

0 


y3/a  Yl-pJ/a 


exp  (-2i/a) 


as) 


d2 

d. 


r  0  P(Y-Z)  >  1 

v  -Z  -  pz2  +  £Y  z1"073  exp  C-2t/P)  P(Y-z)  <  1 

0  0 


(16) 


Rigorous  solutions  to  these  equations  have  not  been  attempted  because 
of  the  extent  of  the  effort  required.  Numerical  methods  must  be  resorted 
to  snd/or  their  properties  may  be  investigated  by  solutions  of  simplex, 
approximate  equation* , 


F.  FLAME  PHOTOMETRY  STUD IKS 

Ths  flams  photo:  itrlc  method  for  the  analysis  of  phoephotux  mi 
of  particular  interest  because  if  the  simplicity  of  use,  ease  of  sample 
preparation,  high  sensitivity,  end  the  possibility  of  coupling  the  flame 
photometer  directly  to  the  Technlcon  for  automatic  analysis  of  the  test 
solutions.  A  preliminary  literature  eaareh  indicated  that  the  fleam 
emission  from  phosphorus  compounds  should  not  bs  a  function  of  molecular 
structure  (Ref.  7).  A  preliminary  experimental  program  was  carried  out 
to  sxamine  the  feasibility  of  the  uithov. 


The  spectra  of  flames  containing  sight  different  orgsnophosphorus 
compounds  were  invastigatad  with  a  Hilgar  Watts  quarts  prism  spectrograph  . 
(Ref.  7).  The  band  heeds  were  the  eame  as  chose  glvan  by  Pearsa  and 
Gaydon  (Ref.  8)  for  PC*  indicating  that  excited  PO  was  formed  by  burning 
all  the  organophosphorous  compounds  invastigatad.  One  of  the  strongeet 
bands  was  in  the  ultraviolet  at  335  ran.  Tha  spectra  appearad  to  be 
structure  dependent,  but  the  qualitatJ/e  nature  of  the  experiment  precluded 
a  definitive  conclusion  at  this  point.  Several  experiments  utilizing  an 
alcholol  burner  containing  tributyl  butyl  phoiphonate  with  a  monochromator  - 
photomultiplier  arrangement  indicated  that  the  relative  intensity  at  335 
run  was  proportional  to  tha  concentration. 
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Th  MV  approach  far  Che  analysis  of  phosphorus  was  discon- 

:  -u'  'I  t- he  basis  01  the  work  of  Brody  and  Chaney  which  indicated  that 
t'-<  it  PO  visible  emission  (526  run)  ed  to  higher  sensitivity  (Ref.  9). 

The  simpie  wick  type  alcohol  burner  was  not  useful  for  quanti- 
.  cv.  Btuiies  at  low  concentrations  so  a  commercial  Jarrel  Ash  burner 
was  purchased.  The  Jarrel-Ash  "Hetco"  burner  Is  an  external  mix  total 
consumption  burner  designed  on  lntr aeon* 1  gnj  flow  principle.  The  ve¬ 
locity  of  the  oxidizing  gas  (air)  passing  the  capillary  orifice  produces 
a  vacuum  in  the  capillary.  The  vacuum  is  used  to  draw  a  sample  up  the 
capillary  aepiracor  and  into  the  shear  stream  of  oxidizing  agent.  The 
detector  consists  of  a  multiplier  phototube  (Sll  response)  and  r  inter¬ 
ference  filter  with  a  peak  transmission  at  526GA  (11 Oa  bandwidth  at  half 
maximum).  A  potential  of  950  volts  was  applied  to  tbc  PM  and  the  algnal 
was  amplified  vi'h  a  Hewlett  Packard  Modal  425  microvolt  meter. 

T!  purpose  of  thic  investigation  was  to  determine  the  lower 
limit  of  detection  of  phosphorus  in  organophosphorous  compounds .  Methanol 
and  isopropanol  were  used  as  '•olventa  because  of  their  low  background 
emission  at  5260.®,.  When  the  flow  rauzf;  of  hydrogen  and  air  were  adjusted 
to  give  the  maximum  signal  in  either  ‘’dvent,  the  calibration  curves  were 
identic*!.  The  Iowa r  limit,  cvf-  of-  phosphorus.'  (30  ppm  solvent) 

(Figure  11)  was 'jecruaied  ,  eb  0,5  pp,  by  buejsfug  out  /fhe  aOlvisnt  Mia o ion 
and  increasing  che  aspiration  rara  of  tha  solvent,  l'hls  is  considerably 
better  than  tha  results  .  -ch  caft  be  obtained  with  emlsaion  spectroscopy 
(50-100  ppm)  and'  ompsrable  to  conventional  flame  photometry  results  of 
1  ppm  of  phi  (phorua  ( inorganic-''- 

The  lower  limit  of  detection  might  be  fnc'Tsaasd  further  by 
using  a  concsve  mirror  behind  the  flame  with  its  centsr  of  curvature  in 
th*  flame  and  by  placing  a  lens,  system  between  the  filter  and  the  PM 
tuba  to  Ln^.raase  th*1;  colli  it  ion  efficiency. 

Thf  amission  spectrum  oi  Bis  in  a  hydrogen  rich  flame  vaa 
■taken  with  a  Perkin-Elmef  Model  99C  scanning  spectrometer  (0.5  an  slits). 
The  dacuctor  waa  an  SMI  935FQ  photcssuJ  tipliac.  The  apoctrum  (Figure  12) 
uncorractsd  fir  spectrometer  respchie,  shows  that  the  band  spectra  of 
die  is  nut  well  resolved.  This  say  be  attributed  to  the  presence  of  the 
solvent  and  the  vise  problem. 

This  vizJble  flame  photometric  method  for  phosphorus  doas  not 
appear  to  be  satisfactory  ">r  total  phosphorus  analysis  because  the  526  nm 
emission  Is  a  function  of  moleculnr  structure.  This  was  first  noticed 
after  the  UV  epcctra  taken  witn  the  HiJser  Watts  Spectrograph  were  care¬ 
fully  ar.aly  ed,  At  similar  phosphoruj  concc ntritione  and  approximately 
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ORAM  OF  Btl/ml  OF  ISOPROPYL  ALCOHOL 


Calibration  curve  for  Bis  using  flame  photcoetric 
detector. 


spertrtm  of  Bis  in  isopropyl  alcohol 


the  same  conditions,  the  darkening  of  the  fiim  varied  with  different 
orgonophos  horous  compounds.  Although  these  data  arc-  only  qualitative, 
they  indicate  chat  the  amount  of  PO  produced  is  influenced  hy  structure 
Brody  and  Chaney  found  a  variation  of  emission  at  526  tun  for  two  poe*i- 
cides  of  different  structure.  At  526  rtm  the  amission  is  due  to  the 
following  preces;  (Rof,  10). 


u  ^  pn  |  u 


HPC  +  hv 


Since  the  amount  of  PO  formed  varies  with  structure,  the  emission  in  a 
hydrogen  rich  flame  at  526  ran  would  also  vary.  This  method,  although 
very  sensitive  and  convenient  will  not  he  accepted  for  the  assessment 
system  because  of  the  dependence  on  structure. 

G.  HASS  DECAY  INVERSION  ANALYSIS 

Tiie  opprcjchus  that  have  been  ucad  to  deduce  the  initial 
particle  •lii'e  d  is  ti  iliti  t  ion  from  measurement  of  mass  decay  (grnv! rational 
settling)  in  a  stirred  chamber  can  generally  be  grouped  into  two  clashes.. 

The  first  nasF  is  an  empirical  approach  and. is  characterized 
by  the  fret  that  cl,  form  of  the  particle  size  distribution  is  assumed 
a  prior  I  and  then  tested  Cor  goodnaa*  of  fit  to  the  observed  mass  decay 
numerical  method*  ua i ng  a  suitable  minimization  criterion,  This  type  cf  . 
solution,  provides  ou.mtitative  cetimate  of  the  gross  pnrticle  statistics 
plus  <•.  qtuillu:  l\*o  picture  of  the  initial  mass  distribution  which  char¬ 
acterizes  the  observed  decoy  measurements . 

The  empirical  approach  locks  the  property,  of  ur.lqueaej* ;  l.e., 

sever,' l  distributions  may  give  approximately  equal  good  fits  to  the 
observed  data.  Hence  ,  the  conclusion  cannot  bs  drawn  from  these  results 
that  the  assumed  distribution  which  yields  the  bast  fit  is  the  general 
solution.  Experience  has  shown,  nowevor ,  that  assumed  mass  distributions 
which  fit  squally  well,  in  general,  will  exhibit  similar  gross  features 
in  the  solution  even  though  derived  from  differing  mathematical  forms. 
Thus,  the  empirical  approach  provides  a  method  for  determining  these 
features  and  tho  choice  of  the  assumed  distribution  becomes  a  matter  of 
convenience. 

A  second  approach  to  the  analysis  of  mass  decay  data  dots  not 
involve  an  assumption  on  the  form  of  the  initial  particle  distribution 
but  attempts  to  derive  it  directly  by  solution  of  the  integral  equation 
describing  the  mass  decay  process,  Ths  solution  is  developed  nurnerically 
by  utilizing  an  appropriate  quadrature  formula  to  obtain  a  set  of 


n imultaneous  linear  equations  which  then  may  be  solved  by  standard  matrix 
inversion  techniques.  The  difficulty  with  this  app.uach  is  that  the 
matrix  to  be  inverted  Is  generally  near  singular ,  and  even  small  errors 
in  the  data  will  produce  unreliable  results.  To  overcome  this  difficulty, 
two  modification#  to  the  basic  approach  have  beon  attempted;  one  is  to 
apply  «  smoothing  operator  to  reduce  the  effect  of  errors  end  the  other 
is  to  oeveiop  a  nonimesr  solution,  loth  or  those  techniques  will  be 
discussed  in  jjrca ter  detail  in  the  following  sections. 

1 .  Empirical  Analysis 

The  lose  of  material  as  a  function  of  time  in  a  chamber 
due  to  gravitational  precipitation  under  conditions  of  continuous  stir¬ 
ring  esn  be  expressed  as 


-  /  exp  r-Cd2T]  f(d)  bd 

a 

where 

M(t) 

Si 

mass  concentration  at  tine 

M(0) 

maul  concentration-  at  tiaa 

T  “  0 

d 

■ 

particle  diameter 

f(d) 
exp  [-Cd2:J 

m 

mesa  distribution  function 

9 

decay  function  for  Stokes' 

settling 

C 

M 

e  constant 

fnd  s  ind  b  define  a  range  of  expected  disasters.  It  remains  to  specify 
the  form  of  f(d)  end  to  e^Juet  the  pereaetere  of  the  aaeuaad  distribution 
to  fit  the  experimental  data. 

Previous  studies  have  assumed  that  f(d)  can  be  represented 
by  the  log-noroel  distribution.  The  work  conducted  on  empirical  methods 
during  the  current  reporting  period  consisted  of  testing  several  other 
distribution  functions  and  comparing  these  with  the  solutions  obtained 
for  the  log-normal. 

One  of  the  functions  testad  la  the  Rayleigh  (Ref.  19)  probability 
distribution.  This  distribution  la  e  special  case  of  the  bivariate 
normal  distribution  and  is  frequsntly  employed  in  radio  caauunlcatlons 
engineering .  The  distribution  is  unimodal  end  positively  skewed;  e 
single  parameter  characterises  the  curve.  The  principal  advantage  of 
this  function  to  the  current  application  it  that  its  form  permits  direct 
integration  of  Equation  (17).  The  single  parameter  then  can  be  estimated 
from  the  observed  data  points. 
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A  second  model  tested  was  based  on  the  gum*  probability 
distribution  function.  This  Is  a  tvc-pareneter  distribution  which  is 
unimoUal  and  can  assume  both  symnetrlc  and  skew  forms  depending  on  the 
^alue  of  the  parameters.  The  exponential  function  is  a  special  case  of 
the  gaum*  distribution.  Using  the  gamma  distribution  in  Equation  (17) 
Coes  not  permit  direct  integration.  neceasiiai.l.ig  m  tiisl  arid  i.»C» 
method  for  estimating  tha  bast  fit  paraoater  values. 


Trees  conducted  with  theta  models  at  wall  at  those  with 
the  log-normal  distribution  indicated  that  approximations  to  s  solution 
could  bs  dsrlvsd  but  that  thsss  gsnarally  satisfied  one  segment  of  the 
observed  data  better  chan  another.  This  is  a  consequence  of  the  char¬ 
acteristic  of  the  test  data  which  indicates  a  rapid  depletion  of  large 
particles  from  the  aerosol  cloud  during  the  initial  period  after  die- 
semination  followed  by  a  residual  decay  of  tha  smaller  particles.  Uni- 
modal  models  for  the  initial  mass  distribution  including  those  discussed 
above  do  not  possess  sufficient  flexibility  to  satisfy  both  of  these 
features  simultaneously.  As  a  result,  further  research  with  empirical 
methods  wag  directed  at  developing  blmodal  mass  distributions. 


The  use  of  bimodal  distributions  for  studying  coagulating 
aerosols  has  bean  reported  earlier  by  Della  Valla,  Qrr  and  Slacker  (Ref. 
7).  The  analytic  form  of  this  modal  w«s  not  coovsniant  for  tha  currant 
application  and  a  stapler  function  was  desired.  For  this  pvrpoar,  *  ' ' 
three-parameter  bimodal  model  was  synthesised  from  the  sum  of  two  Rayleigh 
distributions.  This  sum  can  ba  expressed  as 


f  (d) 


exp  C-dW]  +  (l.K)  4-  axp.Ud (18) 


where  and  w.  era  parameters  of  tha  distributions ,  K  and  (1-K)  are  the 
fractional  mix  of  each.  Hanes,  there  ere  three  parameters  to  be  epaolfVed 
namely  w,  ,  ui.  and  K.  Substituting  in  Equation  (17)  for  f(d)  and  integraeing 
results  In 


MUL  - 


M(0) 


K  exp  [-d2  Ct  +  D/guj] 

2oj2  Ct  +  1 

(1-K)  sxp  [-d2  (2u?  Ct  +  l)/2tt>2] 

- j— — k - 

2^2  Ct  +  1 


(19) 


Assuming  that  d.  is  sufficiently  large  so  that  the  value  of  this  function 
is  negligibly  small  at  the  upper  limit  and  that  is  sufficiently  small 
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»o  that  the  exponential  terra*  approach  unity  at  the  lover  limit,  the 
final  result  becomes 


nvu; 


-  .  *  -  +  _ Q-iQ..,, 

m  '  r. 

2w‘  Ct  +  1  2uj‘  Ct  +  1 


(20) 


This  last  assumption  may  require  values  of  d,  and/or  d_  which  are  soma* 
what  outside  the  Stokes'  Law  region,  but  tha  resultant  error  will  be 
neglected  for  the  tiae  being. 


This  model  was  tested  on  the  same  observed  mass  decay 
date  previously  used  to  test  the  unimodal  models.  Two  algorithms  were 
used  to  estimate  the  model  parameters  in  Equation  (20).  One  of  these 
is  the  classical  Nevton-Rsphson  method  for  the  solution  of  nonlinear 
equations  (Rtf.  12).  T or  ths  currsnt  application,  this  techniqua  was 
constrained  to  seek  values  for  ths  model  perimeters  which  minimised  the 
root-mean-equare  (rms)  error  between  the  obierved  data  and  the  model  values. 
Similarly,  a  more  recent  iteration  technique  developed  by  Marquardt  (Ref. 

13)  was  used,  subject  to  the  same  minimization  criterion.  Both  of  these 
msthoda  vers  programmed  for  automatic  calculation  by  the  1M  1620  computer. 
The  availability  of  both  s.chtaes  provided  en  alternative  when  ecsble 
solutions  were  difficult  to  obtain.  They  alio  provided  a  means  for  vali¬ 
dating  that  the  given  solution  represented  a  global  nna  minimum  and  not 
just  a  local  minimum. 


The  biaodal  model  was  used,  among  othar  things,  to  compart 
the  results  of  «  test  series  using  both  Bis  and  VX  fill  with  a  No.  22 
pltitlc  sphere  as  tha  disseminator.  Ths  rtsultt  from  those  inversion 
snslycas  for  tha  mass  distribution  in  ths  initial  aerosol  cloud  havs  been 
reported  separately*  and  are  amnsrizad  In  Table  6,  The  data  presented 
give  the  distribution  parameters  for  the  modal,  the  mass  recovery  at  24 
minutes  predicted  by  the  model  and  the  rms  error  of  model  fit  to  the  data. 
The  model  parameter  u>,  indicates  ths  mods  diameter  of  the  small  else 
;  particles  in  the  initial  aerosol  cloud  while  gives  the  mods  for  the 
large  sire  components.  In  savarll  teste  thsse'values  are  given  aa  either 
an  arbitrarily  small  value  (<  10*3)  or  as  an  arbitrarily  large  value 
(>  103).  The  latter  caaa  is  used  to  Indicate  that  the  solution  was  not 
capable  of  discriminating  from  the  decay  measurement*  a  stirred  sattling 
of  the  large  site  component  but  rather  this  materiel  ie  depleted  from  the 
cloud  by  other  loss  processes.  In  tha  former  csss,  the  arbitrarily  email 
else  diameter  indicates  that  ths  model  solution  could  not  detect  any  sedi¬ 
mentation  In  the  email  size  component  over  the  time  duration  of  ths  test. 
The  implication  ie  that  s  very  fine  aaroaol  exists  in  the  chamber  aa  s 


residual  mass. 


it 

Latter  Report  to  Mr.  Donald  Buck,  21  April  1967,  also  included  in  thia 
report  were  selected  reeults  from  tests  employing  other  fills  end/or 
different  dissemination  devices. 
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TABLE  6 

INVERSION  ANALYSIS  RESULTS  FOR  BK  AND  VX  TEST  SERIES. 
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A  comparison  of  che  performance  between  these  two  fills 
shows  (1)  that  the  small  size  component  for  VX  is  larger  on  the  average 
than  that  of  Bia  while  (2)  the  24-ainute  recovery  for  VX  is  smaller  on 
the  average  than  thac  for  Bis.  Both  of  these  effects  have  been  noted  in 
earlier  studies;  the  current  model  results  simply  add  further  evidence 
to  this  effect.  These  results  also  show  a  significant  improvement  in  the 
fit  to  be  observed  data  and  it  was  concluded  thac  the  blsodel  model  seems 
te  provide  a  co'V'enient  and  simple  means  for  analysing  mass  decay  dis¬ 
tributions  which  ere  not  adsqustely  described  in  unlmodel  distributions . 

2.  »S  EgfcJflft&fcU 

-3  3 

Both  Che  10  and  10  eolutlona  discussed  above  probably 
arise  as  a  consequence  of  deficiencies  in  the  date.  In  the  former  case, 
the  scatter  in  the  measurements  may  well  mask  che  slow  sedimentation  of 
a  relatively  small  sics  aerosol  cloud.  This  fact  is  somewhat  evident 
from  ths  data  in  Table  6  which  shows  10*3  eolutions  generally  associated 
with  larger  ms  srrors.  In  che  litter  case,  che  10-3  solution  reflects 
the  difficulty  of  obtaining  reliable  measurements  at  very  short  times 
after  the  explosive  event.  The  seneitivity  of  the  bimodal  model  to 
these  measurement  problems  was  investigated  by  two  sets  of  error  analysis 
calculations . 

The  first  was  an  error  analysis  of  the  inversion  model  to 
test  its  sensitivity  to  the  precision  of  ths  mass  recovery  measurements. 

The  procedure  uasd  a  simp la  progressiva  truncation  of  a  set  of  mass  re¬ 
covery  values  from  *  known  particle  else  distribution.  The  Inversion 
modal  was  then  applied  to  each  sat  of  data  and  the  results  compared  between 
the  derived  distribution  and  ths  known  distribution.  The  exact 
value  for. the  node  of  small  sise  fractions  in  this  test  was  Id  microns, 
the  mode  for  che  large  site  fractions  was  130  microns  and  the  mix  be¬ 
tween  the  two-size  distributions  is  in  the  proportion  35/63  percent. 

These  results  are  presented  below  in  Table  7,  A  comparison  cf  the 
known  values  with  those  in  Table  7  shows  a  fair  dagree  of  stability  aven 
for  two-digit  accuracy.  Also,  the  results  show  that  ths  large  slat 
fraction  so am*  to  b«  mors  sensitive  to  errors  in  the  dete  then  the  other 
two  parameters. 


TABLE  7 

RESULTS  FROM  INVERSION  ANALY8IS  OF  TRUNCATED  DATA. 


No. 

Sisnificsnt  Digits 

“1 

Mods  of  Small 
Slsa  Fraction 

U)2 

Mode  of  Large 
Sise  Fraction 

Mixture  of 
Small  end  Larst 

1 

9.00 

129.03 

32/68 

2 

10,35 

156.64 

36/64 

3 

10.04 

151.69 

35/65 

4 

10.00 

150.01 

35/65 

7 

10.00 

150.00 

35/65 
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Sensitivity  of  Che  solutions  to  early  mass  recovery 
measurements  was  explored  by  analyzing  approximately  i  half  dozen  tests 
wichout  the  assumption  of  100  percent  recovery  at  time  zero,  but  with 
the  alternate  condition  chat  the  mass  decay  curve  satisfy  the  first 
measured  recovery  value.  Generally,  this  first  measured  point  is  obtained 
at  1  minute  after  detonation.  A  comparicon  of  the  results  for  these  two 
initial  conditions  shows  rcasonsble  correspondence  between  the  solutions 
for  chs  small  siee  fraction  but  major  differences  in  speclfieat inn  of  the 
large  particle  site  distribution.  Clearly,  most  of  tha  large  particles 
fall  out  or  are  removed  from  the  cloud  within  the  first  minute,  with 
resulting  lost  of  information  on  their  list  distribution. 

3  .  Automatic  Decs  Plotting 

Toward  the  latter  part  of  the  year,  an  additional  computer 
program  was  developed  to  systematize  the  presentation  of  the  Inversion 
analysis  results.  A  sample  of  the  preliminary  format  is  shown  in  the 
attached  Figure  13.  The  graph  showa  the  observed  mass  recovery  values 
(crosses)  and  the  computed  model  mass  decay  curve  (dots).  The  entire 
figure  was  generated  by  an  off-line  IBM  870  Document  Writing  System  frot. 
program  card  instructions  calculated  by  the  computer.  The  system  requires 
approximately  7  minute*  to  plot  the  graph.  Bsverel  modification*  to 
thi*  initial  format  have  been  suggested  and  will  be  added. 

4.  by.  klMH .  &  !  1)8** 

Tachnique*  for  solution  of  Squat ion  (17)  without  tha  aid 
of  e  prior  aseujr.pt ion  on  tha  anaiytlea.1  fora  of  the  mesa  distribution 
coma  under  the  general  category  of  Integral  equation*  of  tha  first  kind. 

A  brief  bibliography  of  acme  of  the  technique*  vhich  have  bean  presented 
in  the  open  literature  ia  given  in  the  references.  None  of  the**  are 
•ntirely  satiafsetory  for  analysis  of  mesa  decay  data  and  in  one  fora  or 
another  are  plagued  by  jaethematicel  problems.  Ona  of  the  most  coonon  of 
theta  difficulties  is  the  instability  introduced  into  the  solutions  aa  a 
result  of  noiay  data.  Numerical  experimentation  with  linear  inversion 
methods  on  known  distributions  has  shown  that  something  of  the  order  of 
eight-place  accuracy  is  required  in  the  measurements  to  recover  tha  ini¬ 
tial  function.  This  type  of  accuracy  is  rarely  available  in  the  meas¬ 
urement  of  physical  systems  and  certainly  not  to  ba  expsetad  from  mats 
decay  data . 

Inversion  techniques  which  deal  with  problems  of  Instability 
have  been  proposed  among  othara  by  Phillips  (Ref.  IS)  and  Twomey  (kef.  16). 
The  basic  approach  is  to  express  the  Integral  equation  as  quadratures  and 
apply  smoothing  constraints  to  the  inversion  of  tha  resulting  linear  system, 
In  the  cast  of  Fhllllpa,  the  smoothing  constraint  la  on  minimising  the 
errors  in  the  second  differences  of  the  indlcial  function.  Twomey  ex¬ 
pands  on  this  general  approach  by  shoving  how  other  minimization  criteria 
can  be  introduced.  The  general  solution  is  given  by  the  matrix  equation 
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121) 


f 


(A 'A  +  -H) 


(A*M  +  >p) 


w  tie  re 


f  »  solution  vector 

A  ■  2  matrix  of  quadrature  coefficient* 

* 

A  ■  the  transpose  of  A 
'  *  the  smoothing  parameter 

H  ■  a  matrix  determined  by  the  minimization  criterion  (elected 
M  •  the  observed  measurement* 
p  ■  a  trial  solution  if  on*  i*  available 


If  p  ■  0  and  H  is  the  matrix  appropriate  to  the  second  difference  minimi¬ 
zation  criterion,  this  solution  is  equivalent  to  that  given  by  Philllpa 
but  possesses  the  advantages  that  (1)  only  one  matrix  inversion  is 

required  and  (2)  the  A  matrix  la  not  required  to  be  aquare.  The  latter 
of  th*ae  permits  aolution  of  an  overdetarmined  system  where  the  number 
of  observation*  exceeds  the  number  of  unknowns . 


The  utility  of  this  method  wts  tested  by  meens  of  e  set 
of  synthetic  mast  dscsy  data.  These  maaa  decay  data  were  synthesized  by 
assuming  a  log-normal  distribution  for  £(d)  (mass  mean  diameter  -7.5u 
geometric  standard  deviation  •  2.3)  and  calculating  tha  mass  decay  using 
Equation  (17).  The  Inversion  than  wts  applied  to  thssa  calculations  in 
an  attempt  to  rscovev  tha  original  log-noroal  distribution.  Tha  results 
of  this  experiment  art  shown  in  Figure  14  for  various  valuta  of  /.  Cas¬ 
par  Ison  of  tha  Inversion  valutt  to  tha  true  distribution  shows  that  the 
matrix  solution  corresponding  to  7  ■  0.1  has  rscovsred  tht  essential 
features  of  the  original  distribution  but  exhibits  consistently  greater 
orolnst*  values.  For  y  *  1.0  tha  matrix  solution  begins  to  loss  tha 
symmetric  character  c f  the  original  distribution.  It  should  be  noted 
that  selection  of  an  appropriate  value  of  7  is  on  a  trial  and  error  basis 
.nd  was  facilitated  in  tha  currant  axpsrlmant  by  a  priori  knowledge  of 
f(d) •  Ir,  the  general  csss,  selection  of  a  value  Tor  7  becomes  a  mors 
subjective  process. 
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5  ,  Nonlinear  Inversion  Method 

The  question  ot  how  to  treat  nc»sy  data  has  w  !  reiontl/ 
been  treated  in  a  unique  .'uproact.  by  nonlinear  met!  ida  developed  at  GCA 
in  another  *  Fnimrien  (l?  can  be  excreseed  as  t  Laplace  Trans¬ 

form  by  making  rl\e  subs  ti  tut  ions 

u  ■  Cd*  and  f  •  t/t  (22) 

in  which  cs-.e  we  obtain 


where 


Mu]  ■»  /  exp  [-u2]  g(u)  du 
J 


2Co  t 


In  thebe  transformations,  t  la  a  acala  time  charactariatic  for  tha 
experiment  which  can  ba  taken  ae  tha  tiaa  Interval  batman  auceaaaiva 
Measurements .  We  now  assume  tha  function  g(u)  to  b«  compoaad  of  a  finite 
aum  of  Dirac  S-functiona  of  height  AP,  at  the  poaitlona  u  ^ ;  namely, 


Mi,  J  -  )  i  •  0,1 . 2n*l  (23). 

J-l 


where 


Xj  •  axp  [*u^  J 


The  problem  of  deriving  particle  size  distribution*  from  observations 
of  raasr  decay  it  analogoua  in  many  retpecta  to  the  meteorological  problem 
of  inferring  the  vertical  atmospheric  temperature  structure  from  meas¬ 
urements  of  the  upwelling  radiance  intensity.  Tha  similarity  arises  from 
the  fact  chut  in  each  cast  tha  physical  process  can  ba  descrlbad  by  an 
integral  equation  with  an  axpontntial  dacay  function  at  ths  kernel.  The 
meteorological  problem  has  received  considerable  attention  due  to  the  planned 
availability  of  radiance  measurements  from  satellite  sensor  systems. 

King  (Ref.  16)  recently  ha#  diacuaied  a  nonlinear  inversion  technique 
which  addteisei  Itself  to  th*  two-fold  problem  of  uniqueness  and  stability 
of  the  solution  in  the  presence  of  noiey  data.  The  above  section  de¬ 
scribes  the  application  of  this  nonlinear  method  to  the  mas*  decay  problem. 
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!i  .lpprox  imu  t  i on  corresponds  physically  to  the  ns.-.impc  Lon  t  ha  c  the 
ini?  tii.  cloud  is  composed  oi  3  series  oi  monodisperso  aerosols  which 
iii  mixed  1  .:ivcn  proportions.  This  system  of  nonlinear  equations 
arises  .  the  construction  of  Gaussian  quadrature  formulae  and  is  dis- 
ussed  in  tfxi3  on  applied  analysis  (Ref.  17).  The  Ret  can  be  solved 
.niquely  for  the  n  weights  .'P,  and  -lie  n  roots  X<  by  meanr  of  the  Prony 
algor  U tun  an  dtsc-ibcd,  for  example  ,  by  Chandrasennar  (.kci,  lfij  . 

Sinn?  only  positive  real  roots  in  the  0  <  X.  <  1  permit 
physical  interpretation,  all  other  roots  must  be  rejected  ind  attributed 
to  arising  from  noise  ir.  the  data.  In  most  applications  of  this  tech¬ 
nique,  the  weight  assigned  to  these  roots  is  small  and  their  exclusion 
can  be  done  with  reasonable  safety.  The  number  of  roots  which  are 
retained  is  a  qualitative  index  of  the  accuracy  of  the  observations. 

In  a  limited  number  of  cases,  this  nonlinear  technique 
was  used  to  analyze  the  test  data.  The  limitation  on  use  of  this  method 
arises  from  the  requirement  that  che  measurements  must  be  spaced  at  equal 
time  increment*.  Relatively  few  tests  satisfied  thit  criterion  for  suf¬ 
ficient  segment  of  the  experiment  to  wur.ant  use  of  che  technique. 
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X,  PROGRAM  FOR  NEXT  QUARTER 

The  general  compos it  ion  of  the  assessment  system  is  now  firm  and 
laboratory  programs  are  to  be  limited  to  modest  efforts  to, 

anH  pvr.lusrton  of  the  beta  absorption- 

V  “  /  .  -  ‘  • - 

impactor  sampler, 

(b)  Modifications  ol  the  Koyco  light  scalier  lug  particle  counter 
for  its  use  in  the  assessment  prograa. 

(c)  Further  investigations  on  problems  of  particle  charging  with 
respect  to  the  mass  decay  of  contained  aerosols. 

The  major  effort  is  to  be  one  of  system  construction  and  application 
in  the  Allowing  areas: 

(a)  Construction  of  an  automated  sequential  filter  sampling  system 
with  provisions  tor  variaole  s.nnpling  programming. 


(b)  Development  of  a  data  format  compatible  with  computer  evaluation 
oi  the  mass  decay  curve  obtained  from  sequential  f liter  samples.  This 

activity  will  be  aimed  at  mioiaiitng  errors  generated  by  human  inter¬ 
vention  in  the  various  step*  required  to  prepare  the  raw  experimental 
data  for  lti  proceseinj  by  the  computer. 
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